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A comparison of the galaxy peculiar velocity field with the 
PSCz gravity field— A Bayesian hyper-parameter method 
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ABSTRACT 



We constructed a Bayesian hyper-parameter statistical method to quantify the 
difference between predicted velocities derived from the observed galaxy distribution 
in the IRAS-PSCz redshift survey and peculiar velocities measured using different 
distance indicators. In our analysis we find that the model-data comparison becomes 
unreliable beyond 70/i _1 Mpc because of the inadequate sampling by IRAS survey 
of prominent, distant superclusters, like the Shapley Concentration. On the other 
hand, the analysis of the velocity residuals show that the PSCz gravity field provides 
an adequate model to the local, < 70/i _1 Mpc, peculiar velocity field. The hyper- 
parameter combination of ENEAR, SN, A1SN and SFI++ catalogues in the Bayesian 
framework constrains the amplitude of the linear flow to be (3 — 0.53±0.014. For an rms 
density fluctuations in the PSCz galaxy number density erf al = 0.42 ± 0.03, we obtain 
an estimate of the growth rate of density fluctuations fcr%{z ~ 0) = 0.42±0.033, which 
is in excellent agreement with independent estimates based on different techniques. 

Key words: methods: data analysis - methods: statistical - Galaxies: kinematic and 
dynamics - Cosmology: observations - large-scale structure of Universe 



1 INTRODUCTION 

The study of peculiar velocity is a powerful tool to ex- 
plore the large-scale structure of the Universe. In the stan- 
dard ACDM cosmology, gravitational instability causes the 
growth of density perturbations and the emergence of the 
peculiar velocity field. In the regime where the density per- 
turbation is linear, the galaxy peculiar velocity (%), sourced 
by th e underlying density field, can be expressed as l|Peebles! 
Il993h 



-Ho/o 

47T 



d 3 x6 m (x,to) 



(1) 



where Ho = H(to) is the Hubble parameter at the present 
epoch, fo is the present day growth rate (henceforth we drop 
the subscript 0) and 5 m is the perturbation to the underly- 
ing dark matter distribution, i.e. S m — (p — ~p)/~p. Assum- 
ing that the observable galaxy distribution links to the un- 
derlying dark matter distribution through a linear, deter- 
ministic bias factor, <5 g = bS m , one can express the above 
equation by substituting the growth rate of density fluctu- 
ations / with the dimensionless parameter /3 = f/b. Ac- 
cording to Eq. (U| the amplitude of the velocity field scales 



linearly with B and this paramete r fully characterises the 
model velocity field l|Peebleslll993h . The value of /3 can be 
estimated by comparing the model velocities predicted from 
Eq. {TJ to measured peculiar velocities estimated from dis- 
tance indicators. A good match between these two vector 
fields would then constitute an observational test for the 
gravitational instability paradigm and for the ACDM model 
dScoccimarro et all l200ll : iFeldman etafl l200ll : IVerde et al] 
2002). Performing such a test is the aim of this paper. 

Comparisons of v-v require all-sky redshift surveys to 
sample the mass distribution in the local Universe in a 
dense and homogeneous way. For this reason, a large num- 
ber of studies have used the IRAS 1.2 J y and PSCz (Point 
Source Catalogue) r edshift catalogues IjFisher et al. I ll995l : 
ISaunders et aljhood h The latter, which we shall use in this 
work, covers about 85 per cent of the sky, contains 12,275 
galaxies at a mean distance of 7500kms~ 1 , and still rep- 
resents the densest redshift catalogue available to date. In- 
deed, the recent 2MRS (T wo Micron All Sky R edshift Sur- 
vey) K B = 11.75 catalogue (|Huchra et al.|[2012l ) has a larger 
depth and better completeness than PSCz, but its sampling 
within 70 h~ Mpc is sparser. 

IRAS redshift catalogues have been extensively used to 
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perform v-v comparisons to estimate 8. iDavis et alj (1 19961 ) 
compared the gravity field obtained from the 1.2 Jy IRAS 
redshift survey with the Tully-Fisher (TF) peculiar veloc- 
ities of ~ 2900 spiral galaxies in the composite Mark III 
catalogue. The presence of systematic discrepancies between 
the two fields precluded the possi bility of measuring 8. The 
likelihood technique developed bv lWillick et all (|l997l Il998h 
was designed to calibrate separately the different subsam- 
ples that constitute the Mark III catalogue. This allowed 
these authors to eliminate these discrepancies and to esti- 
mate 8 ~ 0.49 from 838 galaxies wit h TF distances. An in - 
dependent comparison performed by Ida Costa et ail (|l998l ) 
using peculiar velocities measured from the J - band TF sur- 
vey o f field spirals fSFI. lGiovanelli et al.lll997l ; lHavnes et al.l 
1999T) found 8 ~ 0.6. The sa me likelihood approach as 



Willick et ail l| 19971 ) was used bv lBranchini et alj i|200ll ) to 



compare 989 SFI galaxies having TF-peculiar velocities with 
the PSCz model velocity field. They found that the a linear 
field with 8 ~ 0.42 p rovides a good match to observations. 
IZaroubi et alj (|2002l ) applied a different technique to esti- 
mate a continuous velocity field from the measured veloc- 
ities of both early and late type galaxies in the ENEAR 
(Ida Costa et all |2000| ; iBernardi et ail |2002| ; IWenger et al.l 
I2003T ) and SFI catalogues. They showed that the PSCz grav- 
ity field matches well this velocity field for 8 ~ 0.51. 

Not all v-v c omparisons use TF peculiar velocities. 
iNusser et ail (|200ll ) compared PSCz with the ENEAR early- 
type galaxies having _ D n -cr velo cities and found 8 ^ 0.5. 
iRadburn-Smith et all (|2004l ) and lTurnbull et all (|2012T ) con- 
sidered peculiar velocities of Type la supernovae. They con- 
sidered two different compilations of objects and found 8 ~ 
0.5 and 8 ~ 0.53, respectively. More recently, 2MASS (Two 
Micron All Sky Survey) galaxies dSkrutskie et alj|200q ) were 



also used to perform v-v comparison. iPike fc Hudson! (|2005l ) 
compared the gravity field predicted from 2MASS photom- 
etry and public redshifts to differ ent peculiar velocit y sur- 
veys and found 8 ~ 0.49. Finally, IDavis et ail (|201l] ) com- 
pared the flow traced by the SFI++ samp le of TF velocities 
i|Masters et al.ll2006l ; ISpringob et alj|2007l ) to that predicted 
by the 2MRS sample of galaxies brighter than K s = 11.25 
and found a best fit value 8 ~ 0.33. 

All these studies assume that 8 is independent of scale. 
Since in the standard gravitational instability framework / 
is a scale- independent quantity, this assumption implies that 
galaxy bias is deterministic and linear. In fact, physical pro- 
cesses related to galaxy formation and evolution result in 
a scale-dependent and stochastic bias on small scales, but 
have little impact on v-v analyses in which the distribution 
of mass tracers is smoothed out on larger scales. 

The scatter among 8 values obtained from different v— 
v comparisons comes from several sources such as: (i) the 
use of different techniques to predict the gravity field and 
compare it with measured velocities; (ii) inadequate mod- 
elling of the flow in high density environment; (iii) the use 
of different mass tracers to model peculiar velocities; (iv) 
the use of different velocity tracers that preferentially sam- 
ple different environments; (v) the possible systematic er- 
rors in the calibration of the distance indicator; (vi) pos- 
sible systematic errors i n the use of the distance indicators 
like t h e Malmquist b i as llMalmauisd 1920 ;ILvnden-Bell et all 
1 19881 ; iHudsonl Il993 ; IStrauss and Willickl 1 19951 ); and (vii) 
systematic biases in the model gravity field, e.g. the so- 



called 'Kaiser rocket e ffect' i|Kaiserlll989l ; lNusser et al.ll201ll ; 
iBranchini et al.ll2012l ). 

These considerations together with the improved qual- 
ity of the peculiar velocity data, better understanding of 
the systematics and advances in the modelling techniques, 
lead us to re-examine the issue. In this paper we make an 
extensive comparison between peculiar velocity data from 
different distance ind icators and the PSCz density field 
(Branchini et al.l 1 1999T ) . The rationale behind using differ- 
ent peculiar velocity catalogues is manifold. First, by per- 
forming independent v-v comparisons we can identify and 
correct for systematic errors in the measured velocities. Sec- 
ond, peculiar velocities from different distance indicators 
have different random errors and therefore require differ- 
ent Malmquist bias corrections. Third, possible additional 
sources of systematics can be identified by analyzing the 
residuals in the v-v comparisons restricted to the individual 
catalogues. Finally, we aim at estimating 8 by performing a 
joint comparison that involves different velocity catalogues. 
For this purpose we will use a Bayesian 'hype r-parameter' 
method j|Lahav et al.ll200d ; lHobson et alj|2002l ) designed for 
this purpose. 

This paper is organized as follows. In Section [2] we 
present the PSCz galaxy catalogue and the model velocity 
field observations. In Section[3]we describe the four peculiar 
velocity catalogues considered in this work and the strategy 
we adopt to correct for Malmquist bias. The Bayesian ap- 
proach adopted to perform the v-v comparison is introduced 
in Section 2] and the results, including the estimate of 8 are 
presented in Section [S] Finally, in Section [S] we discuss our 
main conclusions. 



2 MODEL PECULIAR VELOCITIES 

The first ingredient of the v-v comparison is the model 
velocity field. In this work we adopt the one obtained 
i Branchini et all Il999h from the IRAS PSCz catalogue 
( Saunders et al.M2000l ). Panels (a) and (b) of Fig. [1] show 
the angular distribution of PSCz galaxies in Galactic coor- 
dinates. Panel (a) shows objects within a predicted distance 
of 70/i _1 Mpc. This is the distance within which we shall 
perform v-v comparisons. Panel (b) shows objects in the 
shell covering 70-130 ft -1 Mpc. 

Objects in the PSCz redshift catalogues were used to 
trace the underlying mass density field within 300/i _1 Mpc 
under the assumption of linear and deterministic bias 
l|Radburn-Smith et al.ll2004l ). The model velocity field was 
obtained from the positions of galaxies in redshift space 
according to Eq. (JTJ, using the iterative technique of 
lYahil et alj (|l99ll ). Iterations involve only objects within 
~ 130ft _1 Mpc. At larger distances the sampling becomes 
very sparse and objects are used to model an external grav- 
ity field, but their peculiar velocities are not used directly. 

The iterative procedure is potentially prone to system- 
atic errors. Here is a short description of the main ones and 
how to fix them: 



(i) Incomplete sky coverage. The surface density of ob- 
served galaxies drops abruptly near the Galactic Plane, in 
the so-called Zone of Avoidance (hereafter ZoA), clearly 
seen in panels (a) and (b) of Fig. [TJ Different techniques 
have been proposed to fill this region and their impact 
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Figure 1. (a): full-sky PSCz (d < 70 h~ 1 Mpc) plotted in Galactic coordinates (5263 samples in total). The red (light grey) points are 
moving away from us and the blue (dark grey) ones are moving towards us. The size of the points is proportional to the magnitude of the 
line-of-sight peculiar velocity, (b): same as (a) but for 70 < d < 130/i -1 Mpc (3732 samples in total), (c): same as (a and b) but for the 
SFIH — h catalogue, (d): number distribution as a function of distance for PSCz and SFI++ catalogues smoothed over 



on the final predictions is not large. iHudsonl (| 19941 ) and 
iRadburn-Smith et al.l (|2004 ) estimated that typical induced 
errors on ft are well below ~ 8 per cent, an upper limit 
obtained by filling the ZoA with a homogeneous distribu- 
ti on of fictitious galaxie s. Here we use the filling method 
of lBranchini et alj jl999h that consists of randomly cloning 
the observed galaxy distribution from the nearby areas into 
the ZoA. The results obtained with this method were al- 
most i dentical to those obtain ed using a Fourier-based tech- 
nique (|Branchini et al.lll999l ) and co nsistent with t hose o f 
the Wiener-like filtering analysis of IZaroubi et all {2002). 
Analyses performed using mock PSCz galaxy catalogues 
showed that the filling method used in this work induces a 
spurious bulk flow o f ~ SOkms -1 but has negligible impact 
on the estimated ft (Bra nchini et alj|l999h . 

(ii) Uncertainties in modelling the mass distribution be- 
yond 130ft _1 Mpc. These uncertainties mainly induce a 
dipole-like external field that we remove by computing pe- 
culiar velocities relative to the central observer (i.e. in the 
Local Group frame) and by restricting quantitative analyses 
to objects within 70/i _1 Mpc. 

(iii) Treatment of high density regions. IRAS galaxies are 
preferentially late type and therefore undersample the cores 
of galaxy clusters. In addition, regions around high density 
peaks are 'tripled valued', i.e. the same redshift is observed 
at three different positions along the line of sight to the peak. 
To correct these effects we assign an appropriate statistical 
weight to all galaxies near the location of known nearby 
clusters and use the 'robust procedure ' of lYahil et all l|l99l] ) 
to collapse them to cluster's centres. 

(iv) Linear scaling, v{0) oc ft in Eq. (0. Deviations from 



the simple linear scaling are observed if peculi ar velocities 
are p redicted directly from galaxy redshifts (|Davis et al.l 
l201ll ). However, in the framework of iterative techniques, 
the linear scaling is a good approximations and a llows one to 
reconstruct peculiar velocit ies to high accuracy (|Yahil et al.l 
Il99ll ; IBranchini et al.lll999h . 

(v) Kaiser 'rocket' effect. Since the selection function of 
the catalogue is initially computed in redshift space, the re- 
construction p rocedure is potentially prone to the so-called 
'rocket' effect l|Kaiserl Il989l ; IBranchini et al.ll2012f) . Like in 
the previous case, the iterative nature of the reconstruction 
procedure alleviates the impact of the effect that was cor- 
rected by using mock galaxy catalogues. 

(vi) Non-linear effects. Since the model peculiar veloc- 
ities are reconstructed assuming linear theory, non-linear 
motions need to be filtered out. Non-linear velocities arise 
first on small scales. Effective removal is obtained by fil- 
tering the gravity field on a scale Rj, comparable to the 
mean galaxy-galaxy separation, while modelling the velocity 
fiel d at each step o f iterat ions. The value of Rj determined 
by IBranchini et all (|l999T) is shown in figure 3 of their pa- 
per. After iterations, the model velocity field was further 
smoothed with a top hat filter of radius 5 /i _1 Mpc to obtain 
a uniformly smoothed model velocity field. 

(vii) Scale-dependent bias. In this work we assume that f3 
is scale independent, i.e. that IRAS galaxies trace the under- 
lying mass density field according to a simple liner relation 
on scales la rger than Rj . This a ssum ption has been explicitl y 
checked bv lWillick et alj i|l997l ) and lBranchini et al.l (|200lf ). 
Reducing Rj from 5 to 3/i _1 Mpc, effectively probing scales 
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in which the bias is non-linear, caused a modest (~ 5 per 
cent) increase of j3. 

The final result is a linear model for the peculiar veloc- 
ity field specified at the reconstructed real space position of 
8,995 PSGz galaxies within 130 /i~ 1 Mpc that were not col- 
lapsed into galaxy clusters. Their distribution as a function 
of distance is plotted in panel (d) of Fig. [1] (blue curve). 
The vertical line at 70 ft _1 Mpc indicates the volume consid- 
ered in the v-v comparisons ('Data used' in TableQ} and the 
vertical line at 130 fo -1 Mpc separates objects with predicted 
velocities not used in the quantitative analysis ('Outliers') 
from those used to model the gravitational pull of distant 
structures. 

To compare predicted and observed velocities we need 
to interpolate model velocities at the positions of galaxies 
in the peculiar velocity catalogues. This is done by applying 
a Gaussian kernel of the same radius Rj to the predicted 
3D velocity specified at the position of the PSCz galaxies, 

Vrec(Xj), i.e. 



E^Ll«rec(%')exp 



f smo (-^i) 



EN' 
Li ex P 



(2) 



where the sum runs over the N' PSCz galaxies and Xi is 
the position of the galaxy in the peculiar velocity catalogue. 
Interpolated velocities are projected along the line of sight, 
Wi,smo = Vsmo(xi) ■ fi, in order to be compared to the mea- 
sured ones. 

After correcting for systematic errors, the typical ran- 
dom errors on predicted velocities, estim ated from mock 
PSG z galaxy catalogs, is ~ 130 kms -1 fBranchin i et al.l 
Il999h . These errors are much smaller than those in the mea- 
sured velocities and therefore will be neglected in the v-v 
comparison presented in this work. 



3 OBSERVED PECULIAR VELOCITIES 

Measured peculiar velocities are the other ingredient of the 
v-v comparison. In this work we consider four catalogues 
(ENEAR, SN, SFI++ and A1SN) that we briefly describe 
below. We restrict our attention to these catalogues for two 
reasons. First, they are high quality, recently assembled data 
sets that densely sample the peculiar velocity field in the lo- 
cal Universe. Second, peculiar velocities are estimated from 
different distance indicators. These two features minimise 
the chance of systematic errors and thus allow us to per- 
form a joint analysis, reducing random errors considerably. 

There are add i tional velocity catalogue s (SC, 
Giovanelli et all 1 19981; iDale et alj Il999f>. fSMAC. iHudsonl 



1999al ; lHudson et al.ll2004T l, (EFAR lColless et alj|200l . and 
(Willick[l99j|U that have been recent ly used to estimate th e 



bulk flow in the loca l Universe (e.g. IWatkins et al. I (|2009h ; 
iFeldman etahl (|201Ch ). While these other catalogues can 
be useful to measure averaged quantities like the bulk flow, 
they are too sparse and noisy to improve the constrains 
on /3 from a point-by-point comparison, like in the v-v 
analysis. To verify this prejudice we did include them in our 
analysis and found no improvement in the /3 estimate. For 
this reason they are not included in the analysis presented 
here. The four catalogues we use are described as follows: 





d < 70 


70 < d < 130 


ENEAR 


632 


65 


SN 


72 


27 


SFI++ 


2044 


1187 


A1SN 


126 


104 



Table 1. Peculiar velocity samples. The two columns indicate 
the number of samples within the range d < 70 /i _1 Mpc ('in use') 
and 70 < d < 130/i _1 Mpc ('outliers'). 



(i) ENEAR. This is a survey of Fundam ental Plane dis- 
tances to nearby 697 ear l y-type galaxies (Ida Costa et al.l 
|2000| ; iBernardi et~aTll2002l ; IWenger et afll2003l ). either iso- 
lated or in groups. Typical errors for the isolated objects are 
~ 18 per cent of their distance. The characteristic deptld of 
this sample is 29/i _1 Mpc. 

(ii) SN. This sample consi sts of 103 Type la supernovae 
taken from the compilation of lTonrv et al.l (|2003T ). These ob- 
jects are good standard candles and their distances are esti- 
mated more accurately (~ 8 per cent) than in the previous 
case. The characteristic depth of the survey is 32 /i -1 Mpc. 

(iii) SFlH — K This is the larges t and densest survey of pe- 
culiar velocities available to date (|Springob et alj|2007l ). The 
sample considered here consists of 3456 late-type galaxies 
with peculiar velocities derived from the Tully-Fisher re- 
lation. The majority of these objects are in the field (2675) 
and the rest found in groups (726). Their distribution across 
the sky is remarkably homogeneous, as shown in panel (c) 
of Fig. [T] Their radial distribution (red curve in panel (d) 
of the same figure), looks like a scaled-down version of that 
of the PSGz galaxies. The characteristic depth of SFI++ is 
around 40 /i _1 Mpc. Typical errors are of the order of 23 per 
cent. Unlike for other cases, the estimated peculiar velocities 
in th e catalogue have been already corrected for Malmquist 
bias (|Springob et alj|2007l ). 

(iv) A1SN. This catalogue, also known as the 'First 
Amendm ent' supernovae samp le, contains 245 Type la su- 
pernovae (|Turnbullet al.ll2012T ). It was obtained b y merging 



three data sets: (1) a sample o f 106 objects from Jha et all 
l |2007t ) a nd IHicken et al.1 (120091 ) . (2) a collection of 113 ob- 
jects by IHicken et al.l (|2009t ) , a nd (3) 28 objects fro m the 
'Carnegie Supernovae project' (|Folatelli et al.l [2010). The 
characteristic depth of the whole catalogue is 58ft~ 1 Mpc, 
significantly larger than that of the other catalogues. Typi- 
cal errors are of the order of 7 per cent. 

In Table [T] we list the number of objects is each cata- 
logue. Only those within 70 /i _1 Mpc, indicated in column 2, 
are considered for the v-v comparison. At larger distances 
the sampling of the underlying velocity field becomes very 
sparse and random errors too large. Indeed, when we in- 
clude 'outliers' in the range (70-130) /i _1 Mpc (column 3) 
the scatter among the f3 values obtained from the different 
catalogues increases, hinting at possible systematic errors 
in either predicted or measured peculiar velocities. Finally, 



1 The characteristic depth r of each catalogue is defined as 
the error-weighted depth r = X) n iii„r„/ J] n iu„, where w n = 
1 / (a ?,), <y n is the mea s urement error of line-of-sight veloc- 
ity Watkins et al.l 120091; IFeldman et alj l20ld; iMa et alj 1201 ll; 
iTurnbull et alj|2012t) 
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since we perform the v-v comparison in the Local Group 
frame, velocities in the catalogues that are provided in the 
CMB frame are transformed to the Local Group frame by 
subtracting the line of sight component of the Local Group 
velocity determined from th e CMB dipole v = 6 11 km s -1 
toward (/,&) = (269°, +28°) (|Scott fc SmootJl201Ch . 



3.1 Malmquist bias correction 

Peculiar velocities in these catalogues are potentially prone 
to Malmquist bias (MB) since they were estimated (1) 
through 'forward' application of the distance indicator, 
i.e. by estimating a distance-dependent quantity from a 
distance-independent one and (2) under the assumption that 
the estimated distance is the the best estimate of the true 
one (jStrauss and W illick 1995). Of all catalogues considere d 
here only the SFI++ corrects for MB (|Springob et al.ll2007l) . 
In all other cases we performed our own correction using the 
procedure outlined below: 

Let us consider the probability distribution of true 
distance r given the distance inferred from the dis- 
tance indicator d and its error a |Lvnden-Bell et all 1 19881 ; 
IStrauss and Willicldll995T l: 

r 2 n(r)exp(-Ii^) 
P(r\d) = ^—7 Ht, (3) 

where n(r) is the mass density along the radial direction and 
A = (ln(10)/5)<r — 0.46a is the fractional distance uncer- 
tainty of the distance indicators. This conditional probabil- 
ity function can be used to guess the true distance r of an 
object from its estimated distance d, if the the density field 
along the line of sight to the objects, n(r), is known a priori. 
A popular approach is to assume that n(r) is constant and 
the resulting analytic expression is known as homogeneous 
MB. However, n(r) is far from being constant, even on the 
smoothing scales of the model velocity field, and to correct 
for the homogeneous MB we need to follow a different strat- 
egy- 

The key issue is to model n(r). We do this using the very 
same velocity model of Section O More precisely, we use the 
real-space reconstructed positions of the PSCz galaxies as 
mass tracers to interpolate the mass density field on a cubic 
grid of 192/i _1 Mpc and mesh size 1.5/i _1 Mpc, smoothed 
with a Gaussian filter of 5/i _1 Mpc. The field on the lat- 
tice is then interpolated along the line of sight to each ob- 
ject, galaxies and Type la supernovae alike. The value of 
n(r) along the line of sight is specified at the position of 
21 equally-spaced points, with a binning of 1.5/i -1 Mpc. Fi- 
nally, Eq. ([3]) is used to predict r from d using a MonteCarlo 
rejection procedure. 

The MB correction is applied to all objects in the cat- 
alogues, apart from SFlH — h galaxies. In Fig. [21 we compare 
the measured distance (:r-axis, before MB correction) and 
the estimated true distance (y-axis, after MonteCarlo, MB 
correction). Removing this bias preferentially places galaxies 
at larger distances. The magnitude of the effect is quantified 
by the scatter of the points around the black line. It depends 
on the amplitude of the measured velocity error and there- 
fore it increases with the distance and is smaller for Type la 
SN. The dispersion around the black line is not symmetric 



and indicates that, when averaged over many directions, er- 
rors in the observed velocities preferentially scatter objects 
to larger distances. 



4 THE V-V COMPARISON METHOD 

We are now in the position of comparing observed and model 
peculiar velocities at the estimated 'true' distances to mea- 
sure the value of p. 

We compare the observed peculiar velocities of all ob- 
jects in all velocity catalogues with the theoretical predic- 
tions. The latter were basically obtained from Eq. {TJ and 
therefore are sensitive to the mass distribution traced by 
PSCz galaxies out to 300 fo _1 Mpc. However, since the selec- 
tion function of the catalogue drops beyond 130/i _1 Mpc, 
predicted velocities basically probe the mass distribution 
(and the value of /3) within this range. Observed velocities 
are potentially sensitive to mass inhomogeneities on much 
larger scales. However, the fact that we perform our compar- 
ison in the LG frame essentially eliminates the gravitational 
pull from scales larger than 130/i _1 Mpc. Therefore, the v— 
v comparisons at all points within 70/i -1 Mpc effectively 
probe f3 on scales between the smoothing scale (~ 5 /i _1 Mpc 
Gaussian) and the effective size of the PSOz survey (~ 
130/i~ 1 Mpc). In this regard, v-v comparisons are comple- 
mentary to analyses that are based on mu ltipole decom- 
posit i on of the observed v elocity field (e.g. IWatkins et al.l 
120091 ; iFeldman et alj|20ich which, instead, are sensitive to 
the mass distribution on scales larger than that of the pe- 
culiar velocity survey. 

The error budget in quantitative v-v comparisons is 
dominated by uncertainties in the measured velocity. These 
errors increase with distance and depend on the distant in- 
dicator used and therefore may vary considerably from cat- 
alogue to catalogue. In addition, different data sets are po- 
tentially prone to different systematic errors in the calibra- 
tion and in the application of the distance indicator. If the 
goal is to perform a joint v-v comparison, then these differ- 
ences need to be properly accounted for. This can be done 
by adopting the Bayesian hyper-parameter method, which 
is designed to objectively assess whether different errors in 
different catalogues are properly accounted for. The hyper- 
parameter approach is designed to scale the errors of the 
data sets, and then marginalise over all other parameters 
to obtain an estimate of the relative statistical 'weight' of 
the different data sets. In practice, a conventional \ 2 ls de- 
fined as x 2 = T,i( x ° bs ~ 4 hc (#)) 2 /K 2 ), where x° hs and a t 
are the observed quantity and its measurement error, and 
x\ e is the corresponding theoretical value with parameters 
8. Then the hyper-parameter effectively scales the errors as 
otOi. Therefore, by marginalising over the other free param- 
eters, one can compute the distribution of a, which gives an 
objective diagnostic of whether the data sets are problem- 
atic and hence deserve further st udy of the systematic or 
random errors IjLahav et al,l l2000). 

Joint analyses of different velocity catalogues have been 
recently performed to investigate the bulk flow, or higher 
moments of the cosmic velocity field in th e local Universe 
l|Watkins et al.ll2009l ; IFeldman etal]|2010h . When combin- 
ing various data sets the main issue is the freedom in as- 
signing the relative weights of different measurements. The 
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Figure 2. Inhomogcneous Malmquist bias correction. Left panel: the measured distance d (without MB correction) versus true distance r 
(MB corrected); Right panel: measured velocity (without MB correction) versus true velocity (MB corrected). The n(r) function (Eq. J3]l) 
is interpolated by using PSCz density samples. 



standard way of combining two different data sets (A and B) 
is by minimising th e total \ 2 denned as l)Lahav et al.ll2000l : 
iHobson et al1l2002l ) 

X 2 = Xa + Xb- (4) 

This procedure assumes that one can trust the estimated 
random errors, so that the individual \ 2 statistics have equal 
weights. However, when combining two or more different 
data sets with different errors, one may want to assign dif- 
ferent weights to the individual \ 2 statistics 

X 2 = «Xa + &Xb, (5) 

where a and b are the Lagrangian multipliers that constitute 
the Bayesian hyper-parameters. Therefore, even if the mea- 
surement errors are inaccurate, the hyper-parameters can 
assess the relative weight of different experiments, and hence 
let the experiments objectively determine their own weights. 

In the Bayesian hyper-parameter method framework , 
the a posteriori dis tribution of the parameter 6 is defined as 
<|Lahav et al.ll2000h 

-21nP(0|L>) = X> fc ln X fc, (6) 

where D represents the data, the sum is over all data sets, 
Nk is the number of d ata in each data set and x\ * s the \ 2 
of the ith data se t (see lHobson et aT. | |2002| and Appendix B 
of iMa et al.ll2010l for detailed discussions). 

In this work we define the \ 2 f° r each velocity catalogue 

as 

X 2 (/3,a),=E fi Jj* ) , (T) 

i — 1 ^ % ' 

where a™ ca is the measurement error for the line-of-sight 
peculiar velocity vf 1 ^, and a is the hyper-parameter of the 
catalogue. Model velocities i> smo are normalised to ft = 1 and 
linearly scaled by the free parameter ft according to Eq. Q. 
Then the likelihood function becomes 

C{P,a) ~ oT^e - '**. (8) 

Note that this is the likelihood of each individual cat- 
alogue characterized by its hyper-parameter a that repre- 
sents the (unknown) scaling of measurement error. The dis- 



tribution of ft for the single catalogue can be obtained by 
marginalise the likelihood in Eq. (]8]) over a and vice versa. 

In contrast, Eq. © defines the posterior probability of 
the hyper-parameters and /3 given a combination of different 
data sets D. It is obtained from the combination of likeli- 
hoods of the different catalogues in Eq. ([S]). The distribution 
of ft from the joint analysis is obtained by minimising Eq. (JSJ 
with respect to ft. 



5 RESULTS 

In this section we present the results obtained from the 
hyper-parameter method, investigate the v-v comparisons 
separately for each catalogue and assess the goodness of the 
fit by analysing the correlation among the residuals in the 
comparisons. 



5.1 ft value and hyper-parameters 

The results of the hyper-parameter analysis are summarised 
in the top panels of Fig. [3] Panel (a), on the left, shows 
the posterior probability of the ft value obtained form each 
velocity catalogue after marginalising over the correspond- 
ing hyper-parameter a. Different colours refer to different 
catalogues, indicated by the labels. SFI+- 1- galaxies prefer 
a lower value of ft, whereas the distribution of both Type 
la supernovae catalogues peak at larger values. However, 
the overlap among the different probability distributions is 
significant, indicating that the ft values obtained from the 
different v-v comparisons agree with each other. 

This impression is confirmed in Table [2] where we list 
the values of ft at the peak of the various distributions, to- 
gether with the ±lcr width of their Gaussian distributions 
(column 2). For comparison we also list the ft values ob- 
tained from v-v comparisons based on the same velocity 
catalogues, (column 4 in Table [2}. They are listed only for 
reference, since a quantitative comparison should account 
for the different model velocity fields, comparison techniques 
and objects considered in each analysis. No v-v comparison 
was performed with the SFI++ velocities. 

In Fig. [3] the orange curve that peaks at ft ~ 0.53 is 
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Figure 3. (a): Marginalised posteriori distribution of j3 by applying four different catalogues (d < 70 /i _1 Mpc). The orange line is the 
combined constraint by using the hyper-parameter method, (b): Posteriori distribution for hyper-parameter a, same colour scheme as 
panel (a), (c): Goodness of fit for residual velocity correlation xf- 





(likelihood and 10) 


(3 (residual velocity) 


/3 (reference) 


a (this study) 


ENEAR 
SN 

SFI++ 
A1SN 

Combined (hyper) 


0.54 ±0.022 
0.56 ±0.022 
0.51 ± 0.022 
0.56 ±0.032 
0.53 ±0.014 


0.56 ±0.022 
0.56 ±0.022 
0.54 ±0.014 
0.51 ±0.022 


0.50 ± 0.10 CNusser et al. 2001) 

0.55 ± 0.06 fRadburn-Smith et al. 2004) 

0.53 ± 0.08 CTurnbull et al. 20121 


1.31 ±0.04 

1.74 ±0.14 
i cc+0.03 
1 - oo -0.02 
1 41+0.09 

1 - 4i -o.o8 



Table 2. Constraints on f) and a for different catalogues and combinations. We also list the constraints from other published studies. 



the result of combining all catalogues in the joint hyper- 
parameter analysis. The fact that the distribution largely 
overlaps with those obtained for the individual catalogues 
confirms the consistency among the results and allows us 
to estimate f3 with a ~ 2 per cent error. This constraint is 
remarkably tight, but we have to keep in mind that we are 
assuming here that errors are solely contributed by uncer- 
tainties in measured velocities, and that all other possible 
sources, including cosmic variance, can be safely neglected. 

In panel (b) of Fig. [3] we plot the marginalised distri- 
bution of hyper-parameter a for the four catalogues. Their 
values are listed in column 5 of Table The fact that 
a — Nk/x 2 is larger than unity for all subsets suggests 
that random errors have been systematically underestimated 
by a factor ~ q 1//2 . In our analysis we have assumed that 
the error budget is dominated by uncertainties in measured 
peculiar velocities and have ignored errors in the velocity 



model. We've justified this assumption in Se ction [2] based on 
the mock c atalogue analysis performed by (Branchin i et al.l 
1 19991 , | 2001). Another error source that we have neglected 
is associated to the procedure adopted to correct for the 
Malmquist bias. To check whether this can indeed account 
for the remaining error we have performed 1000 Montecarlo 
realization of the MB correction and evaluate its uncertain- 
ties from the scatter in the value of /3. We have done this 
exercise for the ENEAR catalogue and found that the MB 
correction induces an error ap ~ 0.01 of the same size as 
the one from measured velocities. When the two are added 
in quadrature the total error increases by a factor ~ 1.4 1//2 , 
consistent with the hyper-parameters valued. 

We conclude that errors on /3 are contributed by uncer- 
tainties in the measured velocities and in the MB corrections 
and that, for each velocity catalogue, the latter can be esti- 
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mated scaling the former by the hyper parameter a. All /3 
errors quoted in Table [2] have been estimated in this way. 

5.2 Individual v—v comparisons 

Let us now investigate how well predicted peculiar velocities 
match the observed ones in each catalogue. We do this by 
comparing model predictions to observations on a point-by- 
point basis. The results are the scatter plots shown in the 
left panels of Figs. [3]to[S] One point in each plot indicates 
an object in one catalogue. Observed velocities are on the 
y-axis and predicted velocities normalised to /3 = 1 are indi- 
cated on the x-axis. Error bars represent la uncertainties in 
the measured velocities multiplied by the hyper-parameter 
of the catalogue, as indicated in Tabled Straight lines repre- 
sent the best fit value of f3 from the hyper-parameter analysis 
with a slope f3 also listed in Table(2] Panels on the right show 
the velocity residuals u mcas — v ICC (/3 — 1) * /3best-fit computed 
at the estimated 'true' distance of the object. For those cat- 
alogues that contain a sufficiently large number of objects, 
we break down the comparison by distance and show the 
v-v scatter plots for objects in different spherical shells. 

(i) ENEAR catalogue. In Fig.Uwe show the v-v compar- 
isons for objects in three different distance intervals. In the 
innermost shell (d < 35/i _1 Mpc) the general agreement be- 
tween observed and predicted velocities is quite good. The 
few objects that show significant departures from the best 
fit are typically found beyond 20/i -1 Mpc and have posi- 
tive predicted velocities. Similar r esiduals are a l so see n in 
the ENEAR-PSOz velocity maps of lNusser et al.l (120011 ) and 
are mostly galaxies infalling in the Hydra-Centaurus direc- 
tion and outflowing from the Perseus-Pisces complex. The 
Hydra-Centaurus infall is attributed to the Great Attrac- 
tor, an overdensity originally estimated to have a mass of 
~5x 1O 16 M located at (1, 6, cz) ~ (307°, 7°, 4350 /i _1 Mpc) 
l|Lvnden-Bell et ai1ll988l ). subsequently correc ted to M ~ 
8_x 1O 15 M and (289°, 19°, 3200 fe _1 Mpc) l|Tonrv et all 
2000). The absence of large velocity residuals in this shell 
indicates that this prominent infall is well reproduced by 
the model velocity field. This agreement persists out to 
55/i _1 Mpc, i.e. in the second shell, meaning that the pos- 
sible backside infall to the Great Attractor is also correctly 
predicted by the model. 

Measured peculiar velocities in the outermost shell are, 
on average, larger than model predictions. As anticipated, 
the accuracy of the velocity model decreases with distance 
because distant structures are poorly sampled by galaxies 
in the flux-limited PSCz catalogue. One example is the 
Shapley concentration, a large complex of clusters whose 
dyn amical relevance has been outlined by many authors 
fe.g. IScaramella et al. I I 19891 ; lHudsonlH999bl ; iBranchini et al.l 
1999). However, a poor sampling of this supercluster would 
lead to an underestimate of the predicted velocities, whereas 
here we have the opposite effect. The systematic trend ob- 
served in the scatter plot would rather be explained by 
the poor sampling of low density regions, i.e. large super- 
voids. Alternatively, these negative residuals could reflect 
the fact that our model velocity field relies on linear the- 
ory and therefore tends to overestimate peculiar velocities 
near the peaks of the density field, which in turn at large 
distances could be artificially boosted up by shot noise 



|Branchini et al.ll2000l ). Whatever the reason, it is worth 
pointing out that negative ve locity residuals of the same 
amplitude were also found bv lNusser et al.l (|200ll ) for EN- 
EAR galaxies at similar distances and concentrated in the 
area I ~ 0° and -60° < b < -15°. 

(ii) SN and A1SN catalogues. The v-v comparisons of the 
objects in both catalogues shown in Figs. [5] and [S] indicate 
a good match between model predictions and observed pe- 
culiar velocities. This agreement is particularly impressive 
considering the comparatively smaller errors in the mea- 
sured supernova velocities. This result highlights the fact 
that Type la supernova samples are effective probes of the 
underlying density field and should be considered as a seri- 
ous alternative to galaxy-based, peculiar velocity catalogues. 

(iii) SFI++ catalogues. In this case we show the v-v scat- 
ter plots for objects in five different shells. Errors become 
progressively larger with distance, having little statistical 
significance in the outer shell. 

In the innermost shells we notice the presence of several 
discrepant data points, characterised by extremely large pe- 
culiar velocities, all of them outgoing. No other object in 
any other catalogues or in the more external shells has mea- 
sured velocities of the same magnitude. Note that for SFI++ 
galaxies we did not perform our MB correction, but rely on 
the built-in MB corrections. Therefore, the features of large 
velocities may suggest an inadequate correction for the in- 
homogeneous MB. This is quite challenging since in the very 
local region in which the redshift surveys are used to trace 
the density field are often incomplete. 

Finally, in the outermost shells we see that residuals be- 
come systematically more negative, in analogy to what we 
see in the ENEAR catalogue. Explanations proposed to ac- 
count for that behavior are also valid here. 



5.3 Analysis of the velocity residuals 

The Bayesian, hyper-parameter analysis determines the best 
P value that characterises our model velocity field but can- 
not address the question of whether or not the model pro- 
vide an adequate fit to the observed velocities. Inadequate 
fits typically generate spur ious correlations in the v eloc- 
ity residuals. Here we follow I Willick et alj (|l997l Il998l ) and 
IBranchini et al.l (|200ll ) and look for anomalous spatial cor- 
relation in the v-v residual maps. 

First, we define the normalised velocity residual as 

r f o\ _ ^mcas.i — /3 * V IC c,i{P = 1) , . 

°i\P) — a mea ' V y J 

for each galaxy in the catalogue. This quantity depends on 
the free parameter j3 and, when averaged over all objects 
in the catalogues, is minimised by the corresponding best 
fit values in Table [2] Then, we consider all pairs of galax- 
ies in the catalogue and compute the two-point correlation 
function for the velocity residuals: 

v-m-^E^' (10) 

where N(t) is the number of galaxy pairs within predicted 
separation dij < r ± 5/i _1 Mpc and the sum runs over all 
pairs (i, j). If the model velocity field provides a good match 
with the observed velocities then residuals should be spa- 
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Figure 4. Comparison between measured line-of-sight velocities of the ENEAR catalogue and the PSCz gravity field. Left column: direct 
comparison with f) = 0.54 as the best-fit value (Tabled- Right column: residual velocities, i.e. the reconstructed velocities subtracted 
from the measured velocities. All of the errors here are measurement errors multiplied by the best-fit value of the hyper-parameter, 
a = 1.3 (Table [2j- The three rows correspond to different distance intervals. 
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Figure 5. Same as Fig.2]for the SN catalogue. 
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Figure 6. Same as Fig.|4]for the A1SN catalogue. 



tially uncorrelated. In contrast, correlated residuals indicate 
the presence of systematic effects on some scale. 

Normalised velocity residuals and their correlation func- 
tion have been computed for all objects in the catalogues 
and for values of P in the range [0.1, 1] in steps of A/3 = 0.1. 
As an example in Fig. [!J] we show the correlation function 
of the velocity residuals, tp(r), for all ENEAR galaxies with 
separations d scp < 80ft _1 Mpc. The three panels refer to 
different values of /3 specified by the labels, including the 
best-fit value (central panel). Error bars represent Poisson 
noise, N(t)~ ' 5 . It is clear form the plots that the best fit 
model behaves much better that the two other cases ex- 
plored, in which residuals are significantly correlated (or 
anti-correlated) at almost all separations. Focusing on the 
P — 0.54 case we see that the value of 4>( T ) is consistent 
with the null hypothesis of no correlation over most of the 
distance bins. The excess correlation at small scales reflects 
the fact that the model velocity field is smoothed with a 
Gaussian filter of radius 5 h~ Mpc. The residual correlation 
functions of all other catalogues are similar to those of the 
ENEAR galaxies and therefore are not shown here. 

To quant ify the goodness of fit s we compute the follow- 
ing quantity (jBranchini et al.ll200ll ) 

where £(r) = N(r)ip(r), and Abins is the nu mber of bins in 
which we compute ^fT). IWillick et all l|l997t ) showed that if 
residuals are uncorrelated then £(t) behaves like a Gaussian 
random variable with zero mean and variable N(t), and 
that the x| lS distributed like a \ 2 function with a number 
of degrees of freedom equal to N e g ~ 0.87 x N binB , where 
the factor 0.87, computed from mock catalogues, accounts 
for the correlation among the bins. 

We have computed xf as a function of /3 for all velocity 
catalogues. Results are shown in panel (c) of Fig. [3] Different 
symbols and colours represent different velocity catalogues, 
indicated by the labels. To ease the comparison we set the 
minimum of xf equal to zero for all curves. The \ 2 distribu- 
tions are remarkably similar, with minima very close to each 
other and in agreement with the best fit value of p. This re- 
sult is by no means trivial and indicate that the PSCz linear 
velocity model with P — 0.53 provides an adequate fit to the 
peculiar velocities of very different types of objects from dif- 
ferent catalogues. To quantify the agreement we exploit the 



fact that \\ obeys \ 2 statistics and use this fact to compute 
a formal lcr error from Axf = 1. The results of this exercise 
are listed in column 3 of Tableland confirm the qualitative 
agreement found from the visual inspection. 



6 DISCUSSION AND CONCLUSION 

In this work we used the PSCz galaxy redshift catalogue 
to predict the cosmic velocity field in the local Universe 
within the framework of gravitational instability and assum- 
ing linear theory and linear biasing. The model velocity field, 
which depends on a single parameter /3, is compared to the 
observed velocities of different types of objects in the EN- 
EAR, SN, SFI++ and A1SN catalogues. We restrict our 
comparison to objects within 70/i _1 Mpc, where errors on 
measured velocities are reasonably small and we can trust 
the model velocity field. Great care has been taken in cor- 
recting for the inhomogeneous Malmquist bias. This is done 
by using PSCz to trace the underlying mass density field 
and by adopting a Monte Carlo rejection procedure to sta- 
tistically correct for the MB. Finally, model and observed 
velocities have been compared in a Bayesian framework us- 
ing the hyper-parameter method. This technique is designed 
to estimate the statistical weights of different data set in an 
objective way and allows one to jointly analyse different ve- 
locity catalogues. Here are the main results of our analysis. 

(i) The hyper-parameter v-v comparisons performed for 
each catalogue give /3 values that are consistent within the 
errors. This means that the single velocity model is a best 
fit to all data sets. In addition, the best fit P values for the 
ENEAR, SN and A1SN catalogues agree with those obtained 
from previous v-v comparisons that used the same data sets 
but different velocity models and comparison techniques. 

(ii) The joint v-v comparison performed with the hyper- 
parameter technique using all data sets significantly im- 
proves the accuracy in the estimate of p. The best fit value 
is p = 0.53 ± 0.014. 

This result can be used to set a constraint on the growth 
rate of density fluctuations at z ~ 0, i.e. fo- The value 
of this quantity is sensitive to the expansion history of 
the Universe. Its measurement will allow us to under- 
stand whether the current cosmic acceleration is driven 
by a Dark Energy component or perhaps general relativ- 
ity theory needs modification on cosmological scales. Most 
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Figure 7. Same as Fig. [4] for the SFI++ catalogue. 



of the current estimates have been performed at moder- 
ate redshift by quantifying the apparent anisotropies in 
the clustering of galaxies induced b y redshift space dis 
tortions (RSD. IPercival et all 12003; iTegmark et al J 12006 



Guzzo et al.||2008l;ISong fc Perciyajl^OOgtlBlake et aljboil 
Reid et al.ll2012l ; ISamushia et al.ll2012l ; iToieiro et al.ll2012l ) 



Indeed one of the main goal s of ongoing (|Drinkwater et all 
l20ld ; | Eisenstein et al.ll201l | ; Guzzo et al., in prepara tion) 



and planned ( Schle gel et al 



1201 ll ; iLaureiis et "aLlbOHl ) red 



shift surveys is to measure the growth rate at higher red- 
shift, to increase the z-basel ine for this cosmological tes t. As 
pointed out for example bv lHudson fc Turnbulll (2012), the 
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Figure 9. Correlation function of velocity residuals plotted for = 0.2 (left), 0.54 (middle) and 1.0 (right). The samples plotted are 
those pairs with d Bep < 80 h~ 1 Mpc. 



Comparison fas Reference 

SFI++ vs 2MRS 0.31 ± 0.04 Davis et al. (20111 

SN vs PSCz 0.44 ± 0.06 Radburn-Smith et al. (2004) 

A1SN vs PSCz 0.40 ± 0.07 Turnbull et al. (20121 

6dF (RSD) 0.42 ± 0.05 Beutler et al. (2012) 

Combined vs PSCz 0.42 ± 0.033 This study 

Table 3. Constraints on fag from various ca talogues. In the 
ACDM model, / ~ O.^ 55 <Peebldll97lt iLinderlbOOai . 



additional estimates at z ~ can sharpen the observational 
constraints. Moreover, at z ~ the growth rate can be esti- 
mated using techniques alternative to RSD, hence providing 
an important cross-check among the methods. 

In Table [3] we list the most recent estimate of fas, the 
growth rate normalised to the rms mass density fluctua- 
tion on 8 h~ Mpc scales. The last entry is the value ob- 
tained from our analysis computed as fas = /3o"f al , where 
crf al ~ 0.80 ±0.05 is th e rms fluctuation in the numb er den- 
sity of PSCz galaxies (|Hamilton and Tegma rk 2002). This 
result agrees with ACDM model predictions. All estimates 
in Table [3] are in reasonably good agreement. The largest 
discrepancy is below 2a significance. 

(iii) Hyper-parameters. From the hyper-parameter anal- 
yses of the v-v comparison and after marginalising over /3 
we obtain the hyper-parameters of the different data sets. 
We find that all of them are slightly larger but quite close to 
unity, which indicates that velocity errors are slightly over- 
estimated in all velocity catalogues, possibly reflecting the 
fact that in our analysis we did not include uncertainties in 
the MB correction. 



(iv) v-v scatter plots. The inspection of the v-v compari- 
son on a point-by-point basis from scatter plots reveals that, 
on average, the model velocity field provides a good fit to 
observed peculiar velocities in all catalogues. The agreement 
is especially remarkable for the case of Type la supernovae 
that have much smaller velocity errors and constitutes an 

" important check for the gravitational instability scenario. 
The only exceptions are distant ENEAR and SFI++ 
galaxies, for which model prediction systematically overes- 
timate observed velocities. A similar system atic was seen 
in a previous ENEAR vs. PSCz comparison (|Nusser et al.l 
2001), mainly for galaxies in the region I ~ 0° and —60° < 
b < —15°. Among the possible explanations is the possible 
poor sampling of distant inhomogeneities or inaccurate mod- 
elling of the flow around distant density peaks. Additionally, 
a few nearby SFI++ galaxies have peculiar velocities that 
are very large and cannot be matched by model predictions. 
This might just reflect uncertainties in the original correc- 
tion for inhomogeneous Malmquist bias which, in the local 
volume, is notoriously difficult to model due to the incom- 
pleteness in the parent redshift catalogues. 

(v) Goodness of the fit. To evaluate whether our best fit 
model velocity field is also a good model we searched for 
spurious correlations among velocity residuals. We have con- 
sidered all object pairs with separations up to 80/i _1 Mpc. 
Apart from a positive correlation signal at small separations, 
induced by the Gaussian window used to filter out non-linear 
effects, when we set /3 equal to its best fit value we find that 
the residual correlation function is consistent with zero. This 
result further confirms the adequacy of the velocity model 
and the validity of the gravitational instability picture. 
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Although the study of the peculiar velocity field is 
still very data-limited , future surveys, e.g. the 6dF sur- 
vey (| Jones et al.ll2009|^ . or eventually using the 21 cm line, 
e.g. with the Square Kilometre Arravl may provide rich re- 
sources for the study of large-scale structure and cosmic 
flows. But even with these large data-sets, it will still be 
important to carefully assess the potential systematics. 
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